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ABSTRACT 

We use the method of Press, Rybicki & Hewitt (1992) to search for time lags and time leads between different 
energy bands of the RXTE data for MCG-6-30-15. We tailor our search in order to probe any reverberation 
signatures of the fluorescent iron Ka line that is thought to arise from the inner regions of the black hole accretion 
disk. In essence, an optimal reconstruction algorithm is applied to the continuum band (2-4 keV) light curve 
which smoothes out noise and interpolates across the data gaps. The reconstructed continuum band light curve 
can then be folded through trial transfer functions in an attempt to find lags or leads between the continuum band 
and the iron line band (5-7 keV). We find reduced fractional variability in the line band. The spectral analysis 
of Lee et al. (1999) reveals this to be due to a combination of an apparently constant iron line flux (at least on 
timescales of few x 10 4 s), and flux correlated changes in the photon index. We also find no evidence for iron line 
reverberation and exclude reverberation delays in the range 0.5-50ksec. This extends the conclusions of Lee et al. 
and suggests that the iron line flux remains constant on timescales as short as 0.5 ksec. The large black hole mass 
(> 10 8 M Q ) naively suggested by the constancy of the iron line flux is rejected on other grounds. We suggest that 
the black hole in MCG-6-30-15 has a mass of A/bh ~ 10 6 — 10 7 M Q and that changes in the ionization state of 
the disk may produce the puzzling spectral variability. Finally, it is found that the 8-15 keV band lags the 2-4 keV 
band by 50-100 s. This result is used to place constraints on the size and geometry of the Comptonizing medium 
responsible for the hard X-ray power-law in this AGN. 

Subject headings: galaxies:Seyfert, galaxies:individual:MCG-6-30-15, line:formation, methods:statistical, 
X-ray:galaxies 



1. INTRODUCTION 

The X-rays from active galactic nuclei (AGN) are thought 
to originate from the innermost regions of an accretion disk 
around a central supermassive black hole. Thus, in principle, 
the study of these X-rays should allow one to probe the im- 
mediate environment of the accreting black hole as well as the 
exotic physics, including strong-field general relativity, that op- 
erates in this environment. 

In the past decade X-ray astronomy has begun to fulfill 
that promise. Both EXOSAT and Ginga discovered iron K- 
shell features (including the Ka fluorescent line of cold iron 
at 6.4 keV) in the X-ray spectra of Seyfert galaxies which were 
interpreted as 'reflection' of the primary X-ray continuum by 
cold, optically-thick material in the immediate vicinity of the 
black hole (Guilbert & Rees 1988; Lightman & White 1988; 
Nandra et al. 1989; Nandra, Pounds & Stewart 1990; Mat- 
suoka et al. 1990). It was suggested that this cold reflecting 
material was the putative accretion disk of AGN models. With 
the launch of ASCA and the advent of medium resolution spec- 
troscopy, the iron line in several objects was shown to be broad 
(~ SOOOOkms- 1 FWZI) and skewed (Tanaka et al. 1995; 
Nandra et al. 1997). The overall line profiles are in good agree- 
ment with models for fluorescent line emission from the inner- 
most regions of geometrically-thin black hole accretion disks 
(Fabian et al. 1989). Such data allow us to address issues such 
as the location of the radius of marginal stability, the spin of the 
black hole, and the inclination distribution of various classes of 
AGN (see Reynolds 1999 and references therein for a review of 
these studies). In the current, RXTE era, we can now probe the 
iron line and Compton reflection hump in individual objects in 
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some detail (e.g., MCG-5-23-16, Weaver et al. 1998; MCG-6- 
30-15, Lee et al. 1998, 1999a; NGC 5548, Chiang et al. 1999). 

While these spectral studies have been successful, a com- 
plete picture of the AGN phenomenon is not possible without 
addressing the timing properties. Timing studies are important 
for two intertwined reasons. Firstly, AGN are inherently vari- 
able systems. In general, the variability timescale in a given 
object is seen to shorten as one considers higher frequency ra- 
diation. In the X-ray and 7-ray bands, dramatic variability has 
been seen in many Seyfert galaxies with doubling timescales of 
only a few minutes (e.g. see Reynolds et al. 1995). Although 
it is poorly understood to date, the nature of this violent vari- 
ability is a vital component of any final AGN model. Careful 
characterization of the timing properties, as well as determin- 
ing the observed spectral evolution during dramatic temporal 
events, is required if we are to understand this phenomenon. 

Secondly, timing studies are needed to break certain degen- 
eracies that exist in models which, to date, have only been con- 
strained by purely spectral data. The spin of the black hole in 
MCG-6-30-15 provides an excellent example of such a degen- 
eracy — by fitting the 'very-broad' state (Iwasawa et al. 1996) 
of the iron line in this object with models consisting of a thin, 
disk-hugging corona, Dabrowski et al. (1997) inferred that the 
black hole in this AGN must be almost maximally rotating, with 
a dimensionless spin parameter of a > 0.94. However, by in- 
cluding line emission from within the radius of marginal stabil- 
ity, Reynolds & Begelman (1997) showed that a geometry in 
which the X-ray source is at some height above the disk plane 
can produce the same line profile even if the black hole is com- 
pletely non-rotating. While there are subtle spectral differences 
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between the two scenarios (Young, Ross & Fabian 1998) the 
most obvious way of distinguishing these scenarios is through 
their timing properties. The Reynolds & Begelman (1997) ge- 
ometry predicts substantial time delays between fluctuations in 
the primary power-law continuum and the responding fluctua- 
tions in the iron line. More generally, the reverberation charac- 
teristics of the iron line contain tremendous information on the 
mass and spin of the black hole as well as the geometry of the 
X-ray source (Stella 1990; Reynolds et al. 1999). 

The observational situation is more complex. Lee et al. 
(1999b) and Chiang et al. (1999) have analyzed extensive 
RXTE datasets for MCG-6-30-15 and NGC 5548, respectively, 
in order to study the timing properties and spectral variability. 
In both of these objects, the same pattern of spectral variability 
is seen. Firstly, the X-ray photon index displays flux-correlated 
changes in the sense that the source is softer when it is brighter. 
Secondly, and more surprisingly, the iron line flux was found 
to be constant over the timescales probed by these direct spec- 
tral studies (~ 50 — 500ksec). As discussed by both sets of 
authors, these results are difficult to interpret in the framework 
of standard X-ray reflection models since the breadth of these 
lines indicate that they originate from a small region. It ap- 
pears that some feedback mechanism regulates the amount of 
iron line emission in order to produce approximately constant 
iron line flux. Flux-correlated changes in the ionization state of 
the disk represent one such mechanism (we discuss this in more 
detail in Section 5 of this paper). Unless this feedback mecha- 
nism operates instantaneously, we might still expect variability 
of the iron line flux on short timescales. 

Driven by these motivations, this paper addresses the prob- 
lem of determining causal relationships between light curves in 
different X-ray bands, with particular emphasis on timescales 
shorter than those that can be probed by direct spectroscopy. 
In particular, we use the long RXTE observation of the bright 
Seyfert 1 galaxy MCG-6-30-15 reported by Lee etal. (1999a,b) 
and consider the relationship between the 2^-keV band (here- 
after called the continuum band) and the 5-7 keV band which 
contains most of the iron line photons (and hereafter called the 
line band). An important special case is one in which there is a 
linear transfer function relating one band to the other: 

/CO 
dr^(T)a(t-T), (1) 
-oo 

where a(t) and b(t) are continuum and line band fluxes respec- 
tively, and W is the transfer function. Such relationships be- 
tween bands contain much of the important physical informa- 
tion, such as the reverberation characteristics of the iron line. 

Mathematically, the linear transfer equation can be easily in- 
verted using Fourier methods to obtain, 

*(i) = J- r^e-^^M, (2) 
2tt J_oo b(u)) 

where a(cu) represents the Fourier transform of a(t). However, 
in real situations, a large number of regularly sampled measure- 
ments are required to obtain an accurate deconvolution using 
this simple method. More often, deconvolution is achieved us- 
ing maximum entropy techniques or some other regularization 
method (Home et al. 1991; Krolik et al. 1991). 

Another common approach (and one that is often used with 
less well sampled data) is to compute cross-correlation func- 
tions (CCFs), or some variant thereof which accounts for the fi- 
nite and irregular sampling often encountered in real data. The 



discrete correlation function (DCF; Edelson & Krolik 1988) is 
one example of such a variant. Lee et al. (1999b) apply such 
methods to the observation of MCG-6-30-15 considered in this 
paper and detect both phase and time lags between RXTE bands 
(also see Nowak & Chiang 1999). While these methods are 
powerful, it can be difficult to separate subtle time leads/lags 
from the autocorrelation properties of the data. 

Here, we take an alternative approach which is heavily based 
on the method of Press, Rybicki & Hewitt (1992; hereafter 
PRH92). In essence, we use the correlation properties of the 
continuum band data to reconstruct an optimal continuum light 
curve in which the data gaps have been interpolated. Most im- 
portantly, we also compute the expected deviation of the contin- 
uum flux from the interpolated curve. The reconstructed contin- 
uum band light curve is convolved with a trial transfer function 
and compared with the line band light curve in a x 2 sense. We 
then examine changes in the \ 2 statistic as a function of the 
parameters that define the trial transfer function. 

Section 2 recaps the PRH92 method. This is then applied 
to the RXTE data for MCG-6-30-15 in Section 3. The robust- 
ness and validity of our approach is demonstrated by applying 
this method to simulated data (Section 4). Section 5 draws to- 
gether our results and discusses their implications for the nature 
of this source. In particular, we argue that the black hole in this 
AGN has a mass of only 10 6 — 10 7 M©. In order to explain the 
spectral variability, it is suggested that there are flux correlated 
changes in the ionization state of the surface layers of the ac- 
cretion disk. Section 6 presents a short summary of the results 
and relevant astrophysical implications. 

2. THE PROBLEM AND METHOD OF SOLUTION 
2.1. The optimal reconstruction 

The continuum band light curve is reconstructed from the 
data using the technique of PRH92. For completeness, this sec- 
tion summarizes their method. The reader who is primarily in- 
terested in the application of this method may skip to Section 3. 

Suppose that the true flux of the source at time t is s(t), but 
we measure y(t) — s(t) + n(t), where n(t) is the noise in the 
measurement. In our case, the noise is Poisson in nature. Our 
knowledge of s(t) is further impeded by the fact that the mea- 
surement is only made at a finite number of times t i7 where 
i = 1,...,N. We denote y{U) as yi and refer to this as the 
continuum data vector. 

We seek an optimal reconstruction of s(t) which is continu- 
ous in time, s(t), such that 

{e"(t))^{[m-s(t)] 2 ) (3) 

is minimized for all t. As usual, angle brackets denote the ex- 
pectation value. We impose that s(t) is linear in the data vector 
in the sense that 

N 

i=i 

where qi (t) are a set of inverse response functions that are also 
continuous in time. 

Assuming that the noise is uncorrected with both s(t) and 
itself, PRH92 showed that eqn (3) can be minimized to yield, 

N 

m= (5) 
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Here, 

B l0 = (*(*<)«(«,•)> + (n?)J« (6) 

is the total covariance matrix. To keep the notation concise, 
PRH92 define the correlation statistics: 

Cl = (s(U)s(t)) (7) 
dj ee (8(U)a(tj)) (8) 
C(t) = (s(t)s(t)) (9) 

These functions define what PRH92 call the 'covariance 
model' . The expected variance of the real signal from the opti- 
mal reconstruct in eqn (5) is then given by 

N 

(e 2 (t)) = C(t) - ]T Cl (t)(B-% Cj (t). (10) 

Once the covariance model is known, eqns (5) and (10) de- 
fine the optimal reconstruction of the continuum light curve to- 
gether with a statistic measuring the expected deviation of the 
real signal from the reconstruction. 

2.2. The covariance model 

Here, again, we follow the method of PRH92 to determine 
the covariance model for our continuum data. At this stage, we 
make the assumption that the underlying process is statistically 
stationary so that 

(«(ti)*(*2)> - <s(0)*(t2 - ti)) = C(t 2 - h), (11) 

where C (r) is the autocorrelation function that we have to de- 
termine. This function is related to the first order structure func- 
tion V(t) by 

C(r) - (s 2 ) - V(t), (12) 

and V(t) can be approximated by forming pair- wise estimates 
for all distinct pairs of data points in the continuum light curve, 
and then binning by the time lag of the pairs. We find that the 
analytic form 

fits the structure functions of this paper well. In fact, the recon- 
struction is fairly insensitive to the exact analytic form used to 
approximate the structure function. 

Making the reasonable assumption that C(r) — > (s) 2 as 
r — > oo, our final expression for the autocorrelation function 
is 

C(r) = ( S ) 2 + KnalM - Knal(r). (14) 



3 . APPLICATION TO MCG-6-30- 1 5 

In this Section, we apply the method outlined above to a long 
RXTE observation of the bright Seyfert 1 galaxy MCG-6-30- 
15. 

3.1. The RXTE data 

RXTE observed MCG-6-30-15 for approximately 7 x 10 5 s 
starting on 4-Aug-1997. We retrieved these data from the 
NASA-HEASARC public archive situated at the Goddard Space 
Flight Center. Our data reduction closely parallels that of Lee 



et al. (1999a) who has studied the spectral characteristics of 
this observation. Since, as mentioned in the introduction, we 
are interested in the soft X-ray continuum and the iron line 
band, the Proportional Counter Array (PCA) is the appropri- 
ate instrument for us to consider. Examining the housekeep- 
ing files for this observation reveals that Proportional Counter 
Units (PCUs) 3 and 4 suffer occasional breakdown and shut 
off. Hence, we do not consider data from these units and, in- 
stead, extracted STANDARD-2 data from PCUs 0-2. We applied 
fairly standard faint-source screening criteria to these data: the 
source must be at least 10° above the Earth's limb (ELV>10), 
the source must be located within 0.02° of the nominal pointing 
position (OFFSET<0.02), there must be at least three PCUs on 
(NUM_PCU_ON>2), it has been at least 30 minutes since a pas- 
sage of the South Atlantic Anomaly (TIME_SINCE_SAA>30), 
and the electron background is not too high (ELECTRON0<0. 1). 
After application of these screening criteria, approximately 
3.5 x 10 5 s of 'good' data remain. From these data, 2-4 keV 
and 5-7 keV light curves were extracted using 64 s bins. We 
also extracted the 8-15 keV light curve which we will use in 
Section 3.3. 

The background was estimated using the L7-240 background 
models which are appropriate for faint sources such as AGN. 
Background light curves were computed and subtracted from 
the measured light curves in order to form the final background 
subtracted light curves that we shall use in our study. Figure 1 
shows the continuum band (2-4 keV) light curves that results 
from this procedure. For clarity, the light curve shown in this 
figure has been binned with 256 s bins. 

3.2. Searching for lags and leads 

We now apply the procedure outlined in Section 2 to these 
lightcurves. To begin with, we must estimate the structure func- 
tion for these data. Figure 2a shows a pair-wise estimate of 
the continuum band structure function obtained following the 
method of PRH92. This figure also shows our analytic approx- 
imation which is given by eqns (13) and (14) with 



A = 1.757cps (15) 

r = 2.117 x 10 3 s (16) 

ai = 3.805 (17) 

a 2 = 0.195 (18) 

(s)=3.91cps (19) 



Using this covariance model, the PRH92 reconstruction was ap- 
plied to the continuum light curve using N = 5000 data points. 
A portion of the resulting reconstructed light curve is shown in 
Fig. 2b. 

The next stage in the procedure is to convolve the recon- 
structed continuum band light curve with a trial transfer func- 
tion and compare the result with the line band light curve in 
a x 2 sense. We can then minimize the x 2 statistic in order to 
constrain free parameters in the trial transfer function. We also 
minimize x 2 over multiplicative and additive offsets between 
the continuum and line band light curves, i.e. we set 

/ + OC 
i>(T)a(t -r)dr + K, (20) 
-oo 

and minimize over B and K as well as the parameters describ- 
ing the trial transfer function ip. 
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FIG. 1. — 2^-keV band, 3-PCU, light curve for the 1997-Aug-4 RXTE observation of MCG-6-30-15. For display purposes, a bin size of 256s has been used, 
although 64 s bins are used in the analysis presented in this paper. 
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FlG. 2. — Panel (a) shows the structure function for the 2-4 keV band data from MCG-6-30-15 (squares), together with our analytic approximation (solid line). 
Panel (b) shows a portion of the 2-4 keV light curve (open squares and thin error bars) together with the reconstructed light curve (filled circles and heavy error 
bars). 
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In this work, we choose two trial transfer functions. The first 
represents the case where some fraction / tr of the line band flux 
is a delayed copy of the continuum band with a time delay ttr- 

Mt) = (1 - /tr)<J(*) + /tr<*(t - kr). (21) 

The second represents the case where some fraction / tr of the 
line band flux is a delayed and smeared copy of the continuum 
band flux, where a Gaussian kernel is used: 

Mt) = (1 - /trW*) + ^/^exp (J±L^l) , (22) 

<T tr V 7T V 2(T tr / 

No extrapolations were performed during this procedure. In 
order to avoid extrapolating, the x 2 statistic was calculated us- 
ing a subset of data points. For the trial transfer function ipi, 
only data during times t stalt + t tr ,max < t < t cnd - t tr , ma x 
were used to compute x 2 , where t s tart an d t en d are the times 
of the start and end of the reconstructed continuum light curve. 
For ip2> X 2 i s computed based upon data from times t s tart + 

ttr,max ~t~ 2(Ttr 5 max t <i £ cn d ttr,max 2(Tt r .max- 

Figure 3 shows the \ 2 surfaces and confidence contours once 
this procedure has been performed. When displaying the x 2 
surfaces, we plot log 10 (x 2 - Xmin + 1) m order to highlight 
the topography of the surface near the global minimum in the 
surface. It can be seen that the minimum of the x2 surface cor- 
responds to the two lines / tr = and t t r = 0, i.e. no time de- 
layed component of the line band light curve is detected. Here 
we only show the results for tpi — the ip 2 results are trivial 
(i.e. x 2 surface is completely flat) since the preferred solution 
always have / tr = 0. The best fit values of the multiplicative 
and additive constants are B = 0.78 and K = 0.90. 

3.3. The overall time delays between bands 

By considering the / tr = 1 slice through the x 2 surface pro- 
duced with trial transfer function tpi, we can examine overall 
lags and leads between energy bands. Examining the 2^1 keV 
and 5-7 ke V light curves for MCG-6-30- 1 5 in this way, we find 
that the x 2 slice possesses a minimum at zero lag — i.e. we find 
no evidence for overall time lags or leads between the contin- 
uum and line bands down to 64 s, the bin size of the data. Per- 
forming the same procedure for the 2—4 keV and 8-15 keV light 
curves reveals a one bin offset in the position of the minimum 
(Fig. 4), indicating that the 8-15 keV light curve is delayed by 
~50-100s as compared with the 2-4 keV light curve. 

Lee et al. (1999b) have applied CCF methods to this RXTE 
dataset. By carefully comparing with simulations, they find ev- 
idence that the 7.5-10keV band lags the lower energy bands 
with a phase delay of ~ 0.6. They also find evidence that 
the hard band (10-20keV) lags the softer bands with a time de- 
lay similar to that found in this work. Figure 5 shows the DCF 
for our 2-4 keV and 8-15 keV lightcurves (this is very similar 
to Fig. 17 of Lee et al. 1999b). A small time lag of 50-100 s 
between these two bands is evident. Thus, CCF methods and 
the optimal reconstruction method both suggest a time lag of 
50-100 s between the 2^1 keV and 8-15 keV bands. 

3.4. The meaning of an additive offset 

Our fitting in Section 3.2 clearly reveals the need for an addi- 
tive offset (i.e. non-zero K value) between the line and contin- 
uum band light curves. In other words, the fractional variability 
about the mean level is less in the line band that it is in the 
continuum band. 



The spectral analysis of Lee et al. (1999b) allows this be- 
haviour to be understood in terms of the spectral phenomenol- 
ogy. Firstly, Lee et al. found that on timescales of few x lOksec, 
the iron line flux does not track the continuum flux and, in- 
stead, remains approximately constant. Secondly, it was found 
that there are flux correlated changes in the photon index by as 
much as AT ss 0.2 in the sense that higher flux states are softer. 
Both of these spectral changes will tend to reduce variability in 
the line band as compared with the (softer) continuum band. 

4. APPLICATIONS TO SIMULATIONS 

4. 1 . Constructing the simulated light curves 

In order to assess the significance and robustness of the above 
results, this section describes the application of this method to 
simulations. We tailor our simulation to match the RXTE ob- 
servation MCG-6-30- 15 as much as possible. EXOSAT showed 
that the high frequency fluctuations of MCG-6-30- 15 possess 
a power spectrum of the form / _1 ' 36 . We use this power spec- 
trum with an additional low-frequency cutoff at f c = 10~ 6 Hz: 




We then make a (noiseless) simulated continuum light curve, 
F(t), by summing Fourier components of random phase be- 
tween / m i n = 10~ 7 Hz and / max = 1 Hz, i.e. 

F(t) = / /m " X df P(f) sin[2vr / 1 - </>(/)] (24) 

J /min 

where 4>(f) is a uniformly randomly distributed in the range 
to 2ir for each distinct value of /. 

Without loss of generality, we assume that the line band flux 
possesses the same mean normalization as the continuum band 
light curve. However, in order to mimic the situation found in 
Section 3 as closely as possible, we assume that there is an ad- 
ditive offset between the continuum band and line band light 
curves as well as the convolution a transfer function. In other 
words we compute a (noiseless) line-band light curve using the 
expression: 

/+oo 
dr^ hnzl (r)a(t-T)+K. (25) 
-oo 

Here, ^i^ nz i is the non-zero-lag component of our imposed 
simulated transfer function for which we use a Gaussian: 

*W=^exp(-^) (26) 

where / is the fraction of the continuum flux that is delayed, t 
is the mean time delay, and a is the temporal standard-width of 
the smearing. For concreteness, we set 



A = 0.75 


(27) 


/ = 0.15 


(28) 


K = 0.85 


(29) 


to = 1.0 x 10 4 s 


(30) 


a = 2.0 x 10 3 s. 


(31) 



This value of / is approximately the fraction of the line-band 
flux which originates from the iron line, and hence this simula- 
tion crudely mimics the effect of iron line reverberation with a 
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FIG. 3. — Results for MCG-6-30-15: \ 2 surfaces and confidence contours resulting from applying trial transfer function to the reconstructed continuum light 
curves and comparing with the line band light curve. Surfaces are plotted using logio(x 2 — X^in + •*■) as me ordinate in order to display the topography of the 
region near the minimum. Contours are shown the following levels: \ 2 ~ Xmin = ^ ^-2> ^0> 50, 100, 200. The first three of these contours correspond to 

lcr, 90% and 95% for two interesting parameters and are shown in bold. 
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FIG. 4. — The /tr = 1 slice through the \ 2 surface produced by folding the 2-4 keV light curve through the trial transfer function tpi and comparing with the 
5-15 keV light curve. Note the small (1 bin) delay between the two light curves (with the harder band lagging the softer). 
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FIG. 5 . — Discrete Correlation Function (DCF; Edelson & Krolik 1988) between our 2-A keV and 8-15 keV light curves. Note the asymmetry in the DCF which 
validates our detection of a ~ 50 — 100 s time lag between these two bands using the PRH92 method. 
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FIG. 6. — Panel (a) shows a portion of the simulated light curve described in Section 3 of the text. The simulated data possess 64 s bins, and noise is purely due to 
counting statistics. The squares in panel (b) show the results of applying the PRH92 reconstruction to the simulated data (using N = 3000 simulated data points. 
The solid erratic line shows the 'real' simulated signal. Note how well the reconstruction reproduces the 'real' signal during times where data exists, and brackets 
the signal at other times. 
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10 4 s time delay. Our value of K is set to be similar that found 
for MCG-6-30- 15 above. 

From these 'perfect' noiseless light curves, we formed Pos- 
sion sampled noisy light curves assuming a mean count rate 
of 4 cps in both the continuum and line bands, and using 64 s 
bins. The datagap structure of the real MCG-6-30-15 dataset 
was then imposed on the simulated light curves. 

We now examine our realistic, simulated, data in order to as- 
sess how well we can detect the existence of the imposed lag 
and recover the properties of ^i inz i using our method. 

4.2. Extracting the lag from the simulations 

We use the method of Section 2.1 and 2.2 to form an op- 
timally reconstructed, evenly-sampled continuum lightcurve. 
The covariance model used is given by eqn (13) and (14) with 



,4 = 2.5 cps (32) 

t = 6.23 x 10 3 s (33) 

ai = 0.855 (34) 

a 2 = 1 (35) 

(s) = 3.63 cps (36) 



A total of N — 3000 simulated data points were used to form 
the reconstruction which spans a simulated observation time of 
400 000 s. A portion of the simulated dataset and its reconstruc- 
tion are presented in Fig. 6. Note how well the reconstruction 
algorithm recovers the real signal during the times with data, 
and brackets the real signal during other times. 

Figure 7 presents the x 2 surfaces and confidence contours 
that result from passing the simulated light curves through the 
trial transfer functions tpi and ip2> including minimization over 
any additive offset between the continuum and line band light 
curves. Both trial transfer functions clearly detect the imposed 
lag in so far as a a deep and isolated hole is present in the 
X 2 surface at approximately the right time delay, delay frac- 
tion and delay width. Note that the / tr dimension, which has 
been suppressed in the ip 2 plots, has a value of / tr = 0.16 at 
the global minimum. This demonstrates the power of this tech- 
nique for finding and characterizing subtle time lags or leads 
that are present in such data. 

5. DISCUSSION 

In order to bring structure to the discussion that follows, we 
will summarize the pertinent results from this paper. 

1 . We clearly see reduced fractional variability in the iron 
line band (5-7 keV) as compared with the continuum 
band (2-4 keV). This is the origin of the additive offset, 
K, that was introduced in Section 3. The spectral fitting 
results of Lee et al. (1999b) suggests that this is due to 
a combination of a constant iron line flux and flux corre- 
lated changes in the photon index. 

2. Our analysis finds no evidence for iron line reverbera- 
tion effects. By running a number of simulations, we 
find that any reverberation time delays must be less than 
~ 500 s or greater than ~ 50 ks. Together with the above 
result, this suggests an approximately constant iron line 
flux over these timescales. Thus, we can extend the work 
of Lee et al. (1999b) and infer a constant iron line on 
timescales down to 0.5ksec. 



3. Any overall time lag between the 2-4 keV and 5-7 keV 
band is less than ~ 50 s. However, we do find that the 
8-15 keV band is delayed with respect to the 2-4 keV 
band by 50-100 s. We can use this time delay to obtain a 
rough size scale for the Comptonizing cloud that is pro- 
ducing the hard X-rays. Assuming a coronal temperature 
of ~ 100 keV, it takes approximately 3 inverse Comp- 
ton scatterings for a photon to be boosted between the 
2-4 keV and 8-15 keV bands. Thus, the mean free path 
of a photon is approximately 15-30 light seconds. This 
is a lower limit on the size of the Comptonizing region. 

As we will see, this combination of facts presents problems for 
current models. 

5.1. Simple reflection models 

Initially, let us discuss these results in the light of simple X- 
ray reflection models (e.g. George & Fabian 1991). Assuming 
that variations in the primary flux are not accompanied by gross 
changes in geometry, we expect to observe one of two cases. 
Firstly, if the light crossing time of the fluorescing part of the 
disk is shorter than the timescale being probed by the observa- 
tion, an iron line with constant equivalent width will result (i.e. 
the iron line flux will track the flux of the illuminating primary 
X-ray source). On the other hand, if the light crossing time of 
the fluorescing disk is greater than the timescale being probed, 
a constant flux line will result. 

Within the context of these simple reflection models, we are 
forced to conclude that the light crossing time of the fluoresc- 
ing region is larger than <~ 50 ks. Since the line is relativis- 
tically broad, most of the fluorescence occurs in the central 
r ~ 2QGM/c 2 of the disk. Setting the light crossing time of 
this region to be greater than 50 ks gives a black hole mass of 
Mbh ~ 2 x 10 8 M©. Given such a large black hole mass, the 
accretion rate must be less than 1% of the Eddington rate in or- 
der to produce the observed luminosity of Lboi ~ 10 44 erg s _1 
(Reynolds et al. 1997). Furthermore, the size of the X-ray emit- 
ting blobs must be small, rbiob/?"disk ~ 10~ 3 , in order to pro- 
duce the very small time delays seen between different bands. 
Despite being so small, these blobs must be at large distances 
above and below the accretion disk plane, or else one would still 
see iron line variability as a flaring blob illuminated the patch 
of disk directly beneath. 

To date, there are no dynamical measurements of the black 
hole mass in MCG-6-30-15, and hence such a model does not 
explicitly contradict any data. However, there are several in- 
direct arguments that lead us to reject the inference of a large 
black hole mass in MCG-6-30-15. An independent indicator 
of the black hole mass is possible by estimating the bulge mass 
of the host galaxy and then applying the bulge/hole mass re- 
lationship of Magorrian et al. (1998). The B-band luminosity 
of the SO-galaxy which hosts this Seyfert nucleus is approxi- 
mately me = 13.7 (RC3 catalogue), and this is likely to be 
completely dominated by the bulge since the nucleus is heav- 
ily reddened in the B-band and the galactic disk is very weak. 
Using a Hubble constant of Hq = 65kms~ 1 Mpc -1 , the ab- 
solute B-band magnitude of the bulge is then Mb — 19. Using 
the standard relations (Faber et al. 1997), the bulge mass is then 
Mb u ig e ~ 3 x 10 9 M©. Finally, applying the Magorrian et al. 
(1998) scaling factor between bulge mass and black hole mass 
gives Mbh ~ 1 — 2 x 10 7 M©, an order of magnitude smaller 
than the black hole mass estimate in the previous paragraph. 

There are also X-ray constraints that suggest a black hole 
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FIG. 7. — Results for the simulated light curves: \ 2 surfaces and confidence contours resulting from applying trial transfer functions Vi and V"2 to the re- 
constructed continuum light curves and comparing with the line band light curve (allowing for an additive offset between the bands). Surfaces are plotted 
using logu)(x 2 ~ X^in + as me ordinate in order to display the topography of the region near the minimum. Contours are shown the following levels: 
X 2 — Xmin = 2 -3> 4 -6> 9-2, 20, 50, 100, 200. The first three of these contours correspond to 1 — a, 90% and 95% for two interesting parameters and are shown in 
bold. In both cases, the existence of a deep hole in \ 2 space demonstrates that the imposed lag has been clearly detected and its parameters recovered. 
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mass significantly smaller than 10 8 M Q . MCG-6-30-15 has 
exhibited large amplitude X-ray variability on timescales as 
short as 100 s. However, the dynamical timescale of the ac- 
cretion disk where the bulk of the energy is released is t dyn ~ 
10 5 (M/10 8 M Q )s. Thus, if the black hole really is as mas- 
sive as Mbh ~ 2 x 10 8 M©, large amplitude variability would 
be occurring on timescales as short as 10~ 2 idyn- It is diffi- 
cult to conceive of processes which would give such variability. 
The final X-ray argument against a 2 x 10 8 M Q black hole in 
MCG-6-30-15 comes from the power spectrum derived by Lee 
at al. (1999b) and Nowak & Chiang (1999). By comparing 
the power spectral density (PSD) of MCG-6-30-15 with that of 
NGC 5548 and Cygnus X-l, they estimate that the black hole 
in MCG-6-30-15 has a mass of M BH ~ 10 6 M . 

5.2. More complex scenarios 

Since the application of simple X-ray reflection arguments 
led us to deduce an unacceptably large black hole mass, we 
must examine alternative avenues. Indeed, the spectral fitting 
of Lee et al. (1999b) forces us to consider complications be- 
yond the simple reflection picture — In their spectral fitting, 
they found that the Compton reflection continuum fails to show 
the expected correlation with the iron line equivalent width (in 
fact, they are anti-correlated; Lee et al. 1999b). Very similar 
behaviour is also seen in NGC 5548 (Chiang et al. 1999) 

Ionization of the disk surface is one of the few physical 
phenomenon that can (partially) decouple the strength of the 
Compton reflection continuum from the strength of the iron 
line. Matt, Fabian & Ross (1993) demonstrated that the iron 
emission line is more sensitive to ionization effects than the 
general form of the Compton reflection continuum. In other 
words, patches of the disk with certain (surface) ionization pa- 
rameters can produce a Compton reflection continuum without 
producing appreciable iron fluorescence. 

We use this fact to construct the following simple model. Let 
the X-ray flux illuminating the surface layers of the accretion 
disk be 

F(r) oc F^r- f3 . (37) 

A variety of X-ray source geometries give (5 w 3 at large radii, 
and (3 < 3 as one approaches the innermost parts of the disk. 
Now, the ionization parameter of at the surface of the disk is 
given by 

e-^, OB) 

n(r) 

where n(r) is the density of the surface layers of the disk. We 
suppose that n(r) oc r -7 . Hence, we have 



£ cx F x r 



1-13 



(39) 



Standard disk models (Shakura & Sunyeav 1973) give 7 « 2 at 
large radii, and 7 < 2 near the inner part of the disk. Now, sup- 
pose that there exists a critical ionization parameter £ cr j t above 
which there is no iron line produced. For reasonable values of (3 
and 7, this gives a critical radius r cr j t within which no iron line 
is produced. The total iron line flux expected from the object is 
then given by 



/ 



F(r)dr, 



which is readily manipulated to give 



,i(l-7)/(/3-7) 
X 



(40) 



(41) 



For our canonical values of (3 and 7, this gives i*ii ne cx F^ . 
Thus, this simple model produces an iron line flux which is 
anti-correlated with the flux of the illuminating source. Pro- 
vided a strong Compton reflection continuum can still originate 
from the ionized portions of the disk, this type of picture may 
explain the spectral behavior that we observe. 

One simple prediction of this model is that the velocity width 
of the line profile gets smaller as the continuum flux increases 
(due to an outward migration in the inner radius of the line 
emitting region). Of course, the toy model presented above 
only captures the crudest aspects of the problem. Fully self- 
consistent ionized reflection models must be calculated (taking 
into account the vertical structure of the disk; e.g. see Nayak- 
shin, Kazanas & Kallman 1999) and compared with the data in 
order to test whether the picture sketched here is reasonable or 
not. 

Even if global, flux-correlated changes in the ionization 
of the disk surface are responsible for the observed spectral 
changes, we would still expect reverberation signatures on 
short timescales. We have set upper limits of ~ 500 s on the 
timescale of any reverberation delay. If the black hole mass is 
-Mbh ~ 1 x 10 7 M Q , the light crossing time of the iron line 
producing region is ~ 2000 s, and hence we need to infer a 
disk-hugging corona (with h/r ~ 0.3) in order to be compat- 
ible with the reverberation limits. If, instead, the black hole is 
Mbh ~ 1 x 10 6 M©, the light crossing time of the entire line 
producing region is only 200 s and so the X-ray source geom- 
etry is unconstrained by our reverberation limits. The corre- 
sponding Eddington ratios are <~ 10% and ~ 100% for black 
hole masses of 10 7 M Q and 10 6 M respectively. 

6. CONCLUSIONS 

In this paper, we have used an interpolation method based 
upon that of PRH92 to search for temporal lags and leads be- 
tween the 2^keV, 5-7keV and 8-15keV bands in a long 
RXTE observation of the bright Seyfert 1 galaxy MCG-6-30- 
15. In essence, we use the PRH92 method to compute an opti- 
mal reconstruction of the 2-4 keV light curve in which the data- 
gaps are interpolated across. We then fold this reconstructed 
light curve through trial transfer functions and compare with 
data from the other bands in a x 2 sense. 

Our search for lags and leads was tailored to find reverbera- 
tion effects in the iron line which is thought to originate from 
the innermost regions of the black hole accretion disk. We find 
no evidence for any reverberation, and rule out reverberation 
delays in the range 0.5 — 50ksec. We can extend the conclu- 
sions of Lee et al. ( 1 999b), and infer that the iron line possesses 
a constant flux on timescales on timescales as short as 500 s. We 
also find that the hard band (8-15 keV) is delayed by 50-100 s 
relative to the 2-4 keV band. 

We attempt to put these various results together into a co- 
herent picture for this object. The constancy of the iron line 
flux leads one to consider large black hole masses (in excess of 
10 8 Mq). However, such a large mass is found to be unaccept- 
able from the standpoint of both X-ray variability constraints, 
and constraints based on the mass of the galactic bulge. Indeed, 
using the bulge/hole scaling factor of Magorrian et al. (1998), 
we estimate that the hole has a mass of Mbh ~ 1 — 2 x 10 7 M©. 
Given that this is a more reasonable mass estimate, some mech- 
anism beyond the simple X-ray reflection model must be in- 
voked to explain the temporal variability of the iron line and 
Compton reflection continuum. We suggest that flux correlated 
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changes in the average ionization state of the surface layers of 
the accretion disk may be such a mechanism. While we support 
this suggestion with a toy model, the plausibility if this sug- 
gestion can only be assessed once detailed modeling has been 
performed. 
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